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utions  to  verify  the  accuracy  of  the  numerical  model. 


Thermoelastlc  Stresses  in  Multi-Layered  Media 
in  Non-Unifom  Temperature  Fields 

D.  Salinas,  Y.W.  Kwon  and  M.  Neibert 


CHAPTER  1.  ZMTRODUCTZOM 

Because  of  the  mismatch  in  mechanical  and  thermal  properties, 
thermoelastic  stresses  develop  in  multi-layered  media  subjected  to 
temperature  fields.  In  a  1925  paper,  Timoshenko  [1]  presented  the 
solution  for  bending  of  a  bi-metallic  strip  (a  thermostat)  in  an 
uniform  temperature  field.  Closed  form  expressions  for  bending 
stresses,  in  terms  of  strip  thicknesses.  Young's  modulii,  and 
coefficients  of  thermal  expansion,  are  given.  Timoshenko  noted  that 
the  'distribution  of  shearing  stress  along  the  bearing  surface 
cannot  be  determined  in  an  elementary  way'  and  observed  that  the 
magnitude  of  the  shearing  stresses  may  be  of  the  same  order  as  that 
of  the  bending  stresses,  and  additionally,  normal  stresses  along 
the  interface  exist  near  the  ends  of  the  strip. 

In  a  1944  paper,  Goland  and  Reissner  [2]  determined  the  shear  and 
peeling  stresses  along  the  interface  of  a  cemented  lap  joint 
connecting  two  plates  subjected  to  axial  loads  at  their  ends. 
Temperature  effects  were  not  included.  The  closed  form  analytical 
solutions  show  that  the  peeling  stresses  along  the  interface  and 
the  bending  stresses  near  the  cemented  lap  joint  are  as  much  as 
four  times  greater  than  the  axial  stresses  away  from  the  bonded 
interface.  The  shear  stresses  along  the  bonded  interface  are  but 
slightly  less  than  the  axial  stresses  away  from  the  bonded 
interface.  Peak  peeling  and  shear  stresses  occur  at  the  joint 
edges.  All  stresses  depend  on  the  joint  dimensions  and  the 
stiffness  properties  of  the  plates. 

Subsequent  investigators  considered  generalizations  of  the  problem 
of  multi-layered  media  in  uniform  temperature  fields.  In  a 
comprehensive  article  on  the  subject,  Suhir  [3]  summarizes  the 
works  of  others  as  well  as  his  own  on  the  behavior  of  bi-material 
and  tri-material  assembled  electronic  packages  in  an  uniform 
temperature  field,  or  is  subjected  to  applied  axial  load.  The 
results  show  that  for  the  bi-material  assembly,  the  maximum  peeling 
and  shear  stresses  along  the  bonding  interface  occurs  at  the  outer 
edges.  The  article  covers  fatigue,  fracture,  thermal  shock, 
debonding,  and  thermoplastic  failures.  Equations  for  stresses  are 
obtained  from  several  models.  It  is  shown  that  the  stresses  depend 
upon  the  axial  and  interfacial  compliances  of  the  assembly  layers. 
The  middle  bonding  layer  in  the  tri-material  assembly  may  extend 
over  any  fraction  of  the  upper  and  lower  layers.  The  article 
concludes  with  a  list  of  about  240  references  on  the  subject  is 
provided.  Suhir 's  own  model  is  very  general  and  includes  most 
effects.  However  it  is  restricted  to  the  case  of  multi-layered 
media  in  uniform  temperature  fields. 
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The  present  work  generalizes  the  problem  to  that  of  a  multi-layered 
structure  in  a  transient  and  non  uniform  temperature  field  as  might 
occur  during  the  start  of  an  electronic  device.  The  formulation 
consists  of  two  finite  element  codes,  one  for  determination  of  the 
transient  thermal  field  and  the  other  for  the  determination  of  the 
transient  bending,  shear  and  peeling  stresses  resulting  thereof. 
The  finite  element  stress  formulation  utilizes  a  recently  developed 
element  which  provides  for  axial  and  lateral  displacement 
continuity.  The  results  of  the  stress  code  are  in  good  agreement 
with  existing  solutions  [1,3]  to  problems  with  uniform  temperature 
fields. 
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CHAPTER  2  FORMULATION  OP  THE  MODEL 


The  problem  at  hand  is  to  determine  the  stresses  in  a  tri-material 
body  subjected  to  a  temperature  field.  In  order  to  solve  for  the 
stresses  in  a  tri-material  system,  the  problem  was  partitioned  into 
two  parts,  a  thermal  part  to  determine  a  temperature  field,  and  a 
thermoelastic  part  to  determine  stresses  due  to  a  temperature 
field. 

Formulation  of  the  Thermal  Problem 

The  temperatures  in  the  thermal  part  of  the  problem  depend  upon 
material  properties,  system  geometry,  and  the  convective  boundary 
conditions.  In  particular,  the  properties  which  enter  into  thermal 
behavior  include  the  thermal  coefficients  of  expansion,  and  the 
coefficients  of  thermal  condition.  The  thermal  behavior  is  governed 
by  Laplace's  equation 

V*(JcV4»)  ♦  q'"  *  0  (x,y)6  A 

where  T(x,y)  is  the  temperature  field,  q'"(x,y)  is  the  volumetric 
heat  source.  Equation  (1)  is  supplemented  with  convective  boundary 
conditions  along  the  top,  bottom  and  right  surfaces, 

(2) 

on 


and  the  symmetry  condition 


(3) 


along  the  left  (center  line)  boundary  as  shown  in  Figure  2.1. 
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Figure  2.1  Geometry  of  a  Trimaterial  System 

Equation  (1)  is  transformed  into  a  system  of  linear  algebraic 
equations  using  the  Galerkin  form  of  the  finite  element  method 
(FEM) .  In  accordance  with  the  finite  element  method,  the 
rectangular  domain  of  the  structure  is  partitioned  into  the 
discrete  mesh  shown  in  figure  2.2. 
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Figure  2.2  A  Finite  Element  Mesh  of  the  Trimaterial  Medium 

A  piecewise  linear  approximation  T(x,y)  for  temperature  T(x,y)  is 
constructed  as 


T  ^  T  ^  (4) 

where  T(x,y)  is  the  exact  solution  of  Eq.(l),  f  is  the  FEM 
approximate  solution,  {N}^  is  the  row  vector  of  linear  shape 
functions,  and  {T}  is  the  column  vector  of  nodal  temperatures.  A 
residual  function,  R,  which  is  a  measure  of  the  error  between  T  and 
it's  approximation  f  is  defined  by 

Rix,y)  *  V*(icVD  -  q'"  (5) 

Substituting  Eg. (4)  into  Eg. (5)  yields  the  residual 

R{x,y)  »  V* 3  '~Q"Hx,y)  (*) 

Finally,  the  system  of  linear  algebraic  equations  is  obtained  from 
the  othogonalization  of  the  residual  function  with  each  of  the 
(i  «  l,...,n)  shape  functions.  That  is, 

J^jrRdA  -  [JcV(jr'T)  -  g"'l}dA  «  O  (7) 

where  a  change  from  {  }  notation  to  bold  character  notation  has 
taken  place.  Using  Gauss's  theorem,  the  above  equation  becomes 
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A  -  j‘^VJr[JcV(jr')]dA 


and  the  final  'excitation'  vector  tern,  say  Q,  is 


C  -  J  Bq'>'dA 


(11) 


Thus,  after  combining  the  first  and  third  terms  into  F,  where 

r  •  (B  *  ff)  (12) 

Eg. (1)  has  been  transformed  into  the  system  of  linear  algebraic 
equations 

AT  •  r  <23) 

The  solution  of  these  equations  yields  the  temperatures  at  each 
node  of  the  grid.  These  equations  are  for  the  steady  state  problem. 
For  the  transient  problem,  a  transient  tern  is  added  to  the  right 
hand  side  of  Eq.(l),  that  is, 


V*  ( kVT)  q'"  *  a 

or 


(14) 


where  a  is  the  specific  heat  of  the  material.  Preceding  as  before. 
Eg. (13)  becomes 

CT=-AT*F  '•') 

where  the  C  matrix  resulting  from  the  added  a(3T/3t)  term  .  ^iven 
by 


C  =  juNB^dA 


(16) 


In  this  case,  the  nodal  temperatures  are  functions  of  time,  that  is 
T(t). 


A  brief  description  of  the  FEM  formulation  for  stresses  follow.  In 
Figure  2.3,  each  element  has  six  degrees  of  freedom,  axial 
displacements  at  the  four  corner  points,  and  lateral  displacements 
at  the  two  ends.  An  advantage  of  the  element  is  that  axial  and 
lateral  displacement  continuity  results. 


. 

W,  - ►r- 


Figure  2.3  A  typical  element  with  6  degrees  of  freedom 
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The  axial  displacement  field  u(x,y)  is  assumed  to  be  linear  in  both 
the  axial  and  transverse  directions.  That  is, 

u(x,y)  *  ^ <y>  u/  ■*■^2  “/] 

Superscripts  b  and  t  on  the  nodal  displacements  refer  to  bottom  and 
top  displacements  respectively.  The  linear  shape  functions  and 
are  given  by 

H^(y)  -  ^2<^)  “  ^2*7 

where  1  and  h  are  the  length  and  height  of  a  beam  element.  The 
lateral  displacement  field  is  given  by 
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v(x)  (x)  Vi 

(19) 

The 

strain-displacement  relations  are. 

to'  ''‘f’K-i'.V  -  “‘I 

(20) 

and 

(21) 

Defining  the  axial  displacement  vector  as 

=  (  u*  Ux  Uz  Uz  ) 

(22) 

and  the  stiffness  matrix  as 


tiCa]  *  J^^j‘'{B}E{B}'^dydx 


where  the  B  vector  is, 


and  the  force  vector  due  to  temperature  as. 


(Fg)  s  ^f^^j^{B}EatTdydx 


(23) 


(24) 


(25) 


gives  the  bending  matrix  equations  as, 

[JC,1  (6^  -  (26) 

Equation  (10)  defines  the  bending  behavior.  Behavior  due  to  shear 
is  obtained  as  follows.  Define  the  row  vector  of  displacement 
degrees  of  freedom  as 


{6^’’  =  iui  Vj  U2  Uz  Vj) 


(27) 
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The  shear  stiffness  matrix  is  given  by 


[/:.]  -  J^^J^{B')G{B')dydx  (28) 

where 


is') 


a/fa 


a/fi 


aAfj 

dx  "^~dy 

which  gives  the  equations  for  shear  behavior 

[K,\  «  (Ol 


a/f,  dN, 


(29) 

(30) 


A  reduced  numerical  integration  scheme  is  used  to  compute  the 
stiffness  matrix  due  to  shear.  The  matrix  equations  for  bending  and 
shear  behavior  are  combined  to  give  the  stiffness  equations  for  the 
system. 


In  order  to  understand  the  effects  of  material  properties  and 
geometric  configuration  on  system  behavior,  a  parametric  study  was 
undertaken.  The  study  was  conducted  in  two  parts.  In  the  first 
part,  the  effects  of  material  properties.  Youngs  modulus  (E) ,  and 
thermal  coefficient  of  expansion  (a),  were  investigated.  This  is 
presented  in  Chapter  3.  The  second  part  considered  the  effects  of 
system  dimensions  on  system  behavior.  The  effect  of  system  geometry 
on  system  behavior  consisted  of  two  subparts.  The  first  subpart 
considered  the  effect  of  midlayer  thickness  on  system  behavior. 
This  is  presented  in  Chapter  4.  The  second  subpart  considered  the 
effect  of  midlayer  length  on  system  behavior.  This  is  presented  in 
Chapter  5.  These  studies  were  conducted  for  a  uniform  temperature 
field  over  the  trimaterial  system.  Chapter  6  considers  the  effect 
of  a  nonuniform  temperature  field  over  a  trimaterial  system  with 
application  to  electronic  packages  where  the  midlayer  is  the 
adhesive  that  joins  the  chip  to  a  substrate. 
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CHAPTER  3  EFFECT  OF  MATERIAL  PROPERTIES  ON  SYSTEM  BEHAVIOR 

In  this  study,  the  trimaterial  system  was  subjected  to  a  uniform 
temperature  field  IOC  •>€  above  the  stress  free  temperature.  Each  of 
the  three  layers  had  a  depth  of  1  mm.  The  top  and  bottom  layers  had 
equal  half-lengths  of  10  mm. ,  and  the  midlayer  extened  across  half 
of  this  width  as  shown  in  Figure  2.1  below.  The  system  is  symmetric 
with  respect  to  the  y-axis. 

The  range  of  properties  shown  below  in  Table  3.1  were  selected  from 
the  properties  of  materials  which  comprise  trimaterial  electronic 
packages,  that  is,  ceramic  substrates,  silicon  or  arsenic  chips, 
and  solder  adhesives.  Although  trimaterial  configurations  used  in 
applications  other  than  electronic  packages  may  have  significantly 
different  properties  than  those  considered  here,  it  is  expected 
that  the  observations  made  here  remain  valid.  The  quantitative 
would  be  different,  but  the  behaviorial  patterns  would  not  change. 


VALUE 

E  Youngs  Modulus  GPa 

a  Therm.  Coeff.  ®C 

Minimum 

10 

1x10“^ 

Maximum 

100 

400x10^ 

Table  3 . 1  Range  of  property  values 
Effect  of  Youngs  modulus 

The  first  study  investigated  the  effect  of  Youngs  modulus  on  system 
behavior  when  both  the  temperature  and  coefficients  of  thermal 
expansion  are  uniform  throughout  the  system  and  each  of  the  layers 
have  different  Youngs  modulii.  As  expected,  the  results  show  that 
no  stresses  develop. 


Effect  of  thermal  coefficient 

In  this  study,  the  three  layers  had  the  same  Youngs  modulus,  100 
GPa.  Again,  the  temperature  field  was  uniform  for  all  layers  at  100 
*C  above  the  stress  free  state.  To  determine  the  effects  of  thermal 
coefficient  of  expansion  on  system  behavior,  various  combinations 
of  thermal  coefficients  of  expansion  were  selected  for  the  three 
layers.  The  eleven  cases  shown  in  Table  3.2  below  were  considered. 
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1  Case 

Material  A 

Material  B 

Material  C 

1 

100 

50 

100 

2 

100 

200 

100 

3 

300 

100 

10 

4 

300 

10 

100 

5 

100 

300 

10 

6 

150 

100 

10 

7 

200 

100 

10 

8 

400 

100 

10 

9 

200 

100 

1 

10 

200 

100 

50 

11 

200 

100 

75 

Table  3.2  Cases  for  study  of  effects  of  thermal  coefficient  of 
expansion 

Cases  1  and  2.  In  cases  1  and  2,  outer  layers  A  and  C  have  the 
same  thermal  expansion  coefficients.  For  case  the  midlayer  has 
a  lower  coefficient  of  expansion  than  the  outer  layers,  and  in  case 
2  the  midlayer  has  a  higher  thermal  coefficient  of  expansion  than 
the  surrounding  outer  layers. 

For  case  1,  the  lower  coefficient  of  expansion  of  the  midlayer 
should  cause  the  midlayer  to  expand  along  it's  top  and  bottom 
interfaces.  The  shear  and  normal  stress  results  presented  in  figure 
3.1  shows  that.  These  figures  also  show  that  these  interface 
'stretching'  shear  stresses  are  accompanied  by  interface  peeling 
normal  stresses.  It  is  also  seen  that  the  shearing  stresses,  for 
the  particular  geometric  configuration  previously  described,  are 
three  to  six  times  larger  than  the  normal  stresses.  It  is  also 
noted  that  the  maximum  shearing  and  normal  stresses  occur  at  the 
corner  ends  of  the  interfaces. 

The  results  for  case  2,  wherein  the  midlayer  had  a  higher  thermal 
coefficient  than  the  neighboring  upper  and  lower  layers,  might  have 
been  anticipated  from  the  results  of  case  1.  Whereas  the  midlayer 
was  previously  being  'tensile'  shear  along  its  interfaces,  it  is 
now  subjected  to  the  opposite  effect  of  'compressive'  shear 
stresses  along  its  interfaces.  The  'compressive'  shear  stresses  are 
now  accompanied  by  normal  bearing  stresses  on  the  upper  and  botton 
surfaces  of  the  midlayer.  Again  the  magnitudes  of  the  shearing 
stresses  are  three  to  six  times  larger  than  the  normal  stresses  and 
the  maximum  magnitudes  of  stress  occur  at  the  corners  of  the 
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nidlayer.  In  both  cases  the  naximun  shear  stresses  occur  at  the 
innermost  point  of  the  midlayer  and  the  maximum  normal  stress 
occurs  at  the  outermost  midlayer  point. 

Cases  3.  4  and  5.  These  cases  were  selected  to  assess  the  effect  of 
the  ordering  of  thermal  coefficients  as  follows.  In  case  3, 

>  Oq,  that  is,  the  coefficient  of  thermal  expansion  decrease  from 
top  to  bottom  of  the  trimaterial  system.  In  case  4,  the 
coefficients  from  top  to  mid  to  bottom  layers  were  selected  in  the 
order,  largest  value:  smallest  value:  mid  value.  In  case  5,  the 
order  was  changed  to  mid  value:  largest  value:  smallest  value.  The 
magnitudes  of  the  thermal  coefficients  of  expansion  for  all  three 
cases  are  shown  in  figure  3.2. 

The  results  from  these  three  cases  are  presented  in  the  same 
figure.  A  comparison  of  the  results  shows  that 

•  The  smallest  shear  and  normal  stresses  occur  for  case  3,  where 
the  thermal  coefficients  of  expansion  decrease  from  the  top  to  the 
bottom  of  the  trimaterial  system.  In  this  case,  the  deformation  o 
the  mediiun  is  similar  to  that  of  pure  bending.  As  a  result,  these 
stresses  are  significantly  less  than  the  stresses  obtained  in  cases 
4  and  5. 

•  In  cases  where  the  thermal  coefficient  of  expansion  of  the 
midlayer  is  the  largest  or  the  smallest  compared  to  those  of  the 
outer  layers,  the  interface  shear  and  normal  stresses  become  large. 

•  When  the  thermal  coefficient  of  expansion  of  the  midlayer  is 
the  largest,  the  normal  stress  is  a  bearing  stress  and  the  shear 
stresses  contract  the  midlayer  as  shown  in  case  5.  On  the  other 
hand,  the  smallest  thermal  coefficient  of  expansion  (case  4)  causes 
the  peeling  normal  stress  at  the  interface  and  the  shear  stress  in 
the  opposite  direction  of  that  in  case  5. 

•  Again  in  all  cases,  the  interface  shear  stresses  are 
significantly  larger  than  the  interface  normal  stresses. 

Cases  3.  6.  7  .  and  8.  In  this  study,  the  thermal  coefficients  of 
expansion  of  the  mid  and  bottom  layers  was  held  constant  while  the 
thermal  coefficient  of  expansion  of  the  top  layer  was  varied  as 
shown  in  figures  3.3  and  3.4.  In  all  cases  here,  the  coefficients 
of  expansion  had  values  decreasing  fro  top  to  bottom  layer.  With 
regard  to  the  top  layer  whose  thermal  coefficient  of  expansion 
varied  in  decreasing  order  as  follows:  case  8  (400xl0'V *C) ,  case 
3  (300x10  ®/ ’C) ,  case  7  (200x10  ®/ 'C) ,  and  case  6  (150x10  ®/ ‘C) .  The 
values  of  the  thermal  coefficient  of  expansion  for  the  mid  (B)  and 
bottom  layers  (C) ,  were  fixed  at  100xl0'®/*C,  and  10xl0'®/*C 
respectively. 

The  results  of  this  study  presented  in  figures  3.3  through  3.4  show 
the  following. 
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E^  =  Ep=E(,=  100  GPa 

Figure  3.2  Shear  And  Normal  Stresses  On  The  MicDayer  For  Cases  3,  4,  And  5 
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Figure  3.4  Shear  Stress  At  Interfeces  For  Cases  3, 6,  7,  And  8 


•  The  normal  stresses  are  a  naxinun  when  the  top  layer  has  its 
largest  thermal  coefficient  of  expansion. 

•  These  normal  stresses  decrease  with  decreasing  value  of 
thermal  coefficient  of  expansion. 

•  If  the  thermal  coefficient  of  expansion  of  the  top  layer  is 
much  larger  than  that  of  the  midlayer,  the  interface  normal  stress 
is  a  peeling  stress.  On  the  other  hand,  if  the  thermal  coefficients 
of  expansion  of  top  and  middle  layers  are  not  so  much  different, 
the  interface  normal  stress  becomes  a  bearing  stress. 

•  The  maximum  peeling  stress  occurs  at  the  inner  corner  point  of 
the  upper  interface  surface  and  the  outermost  corner  point  of  the 
lower  interface  surface. 

•  The  nozmal  stress  distribution  does  not  vary  significantly 
away  from  its  peak  corner  values. 

•  The  interface  shear  stresses  are  significantly  larger  than  the 
interface  normal  stresses. 

•  The  shear  stresses  achieve  their  maximum  values  at  the  outer 
corner  for  the  upper  Interface  surface  and  the  inner  corner  for  the 
lower  Interface  surface  except  for  cases  6  and  7.  The  maximiim  shear 
stress  on  the  upper  interface  surface  is  significantly  larger  than 
the  maximum  shear  stress  on  the  lower  interface  surface. 

•  The  shear  stresses  achieve  their  maximum  value  when  the 
thermal  coefficient  of  expansion  of  the  top  layer  is  largest.  The 
magnitudes  of  shear  stresses  decreases  with  decreasing  magnitude  of 
the  coefficient  of  expansion  of  the  top  layer. 

It  should  be  noted  that  the  results  reported  above  are  for  the  case 
where  the  lower  layer  has  a  fixed  value  of  coefficient  of  thermal 
expansion  below  the  values  of  thermal  coefficient  of  expansion  for 
the  mid  and  upper  layers.  If  the  roles  of  upper  and  lower  layers 
were  reversed,  that  is,  if  the  upper  layer  had  a  fixed  thermal 
coefficient  of  expansion  lower  than  the  mid  and  bottom  layers  then 
the  previous  results  would  be  valid  simply  by  substituting  the  word 
'upper'  by  'lower'  and  the  word  'lower'  by  'upper'. 

Cases  7.  9.  10.  and  11.  In  this  group  of  cases  the  effect  of  the 
decreasing  values  of  thermal  coefficient  of  expansion  of  an  outer 
layer  on  system  behavior  was  investigated.  As  shown  in  figures  3.5 
and  3.6,  the  values  of  the  thermal  coefficient  of  expansion  of  the 
top  and  mid  layers  were  fixed  at  200x10'”/ ’C  and  100x10'”/ *C 
respectively,  while  the  coefficient  of  thermal  expansion  of  the 
bottom  layer  was  varied  from  75x10'”/ *C  (case  11),  50x10'”/ *C  (case 
10),  10xl0'”/‘C  (case7)  to  lxl0'”/*C  (case  9).  The  results  of  this 
study,  are  presented  in  Figures  3.5  through  3.6  and  are  briefly 
summarized  here. 
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Figure  3. 5  Normal  Stresses  At  Interfeces  For  Cases  7. 9, 10,  And  1 1 
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Figure  3. 6  Shear  Stress  At  Inter&ces  For  Cases  7, 9, 10,  And  1 1 


•  Both  the  normal  and  shear  stresses  are  significantly  lower 
than  in  the  previous  study  where  the  thermal  coefficient  of 
expansion  for  an  outer  layer  took  on  increasing  values. 

•  Again  the  interface  shear  stresses  are  significantly  larger 
than  the  interface  normal  stresses.  And  again,  the  largest 
interface  normal  and  shear  stresses  occur  at  the  inner  and  outer 
corner  points  of  the  midlayer. 

•  As  the  value  of  the  thermal  coefficient  of  expansion  decreased 
in  the  sequence  case  11,  case  10,  case  7,  case  9,  the  normal 
stresses  at  the  left  side  of  the  upper  interface  surface  were 
always  peeling.  The  magnitude  of  these  upper  left  inner  corner 
peeling  stresses  decreased  as  the  thermal  coefficient  of  expansion 
of  the  bottom  layer  decreased. 

•  For  cases  11  and  10,  the  normal  stresses  along  the  upper 
interface  surface  were  always  peeling  stresses.  However  for  cases 
7  and  9,  the  normal  stresses  continue  to  decrease  along  the 
interface  and  became  bearing  stresses  at  the  outer  corner  of  the 
interface. 

•  Somewhat  similar  behavior  with  regard  to  normal  stresses  was 
observed  for  the  lower  Interface  surface.  Here  all  the  normal 
stresses  were  peeling  at  the  outer  corner,  with  decreasing 
magnitude  as  the  thermal  coefficient  of  expansion  of  the  bottom 
layer  decreased.  At  the  inner  corner  of  the  lower  interface 
surface,  the  normal  stresses  were  bearing  stresses  with  decreasing 
magnitude  as  the  thermal  coefficient  of  expansion  of  the  bottom 
layer  increased. 

•  In  all  cases  the  shear  stresses  along  the  upper  interface 
surface  changed  direction  along  the  interface  from  a  maximum 
'negative'  value  at  the  inner  corner  to  a  maximum  'positive'  value 
at  the  outer  corner.  These  shear  stresses  were  of  a  local  nature 
with  very  sharp  decrease  to  insignificant  stress  magnitude  over 
most  of  the  interface  length.  The  magnitude  of  the  interface  shear 
stresses  decreased  as  the  magnitude  of  the  thermal  coefficient  of 
expansion  decreased. 


•  On  the  lower  interface  surface,  the  shear  stresses  were  of  a 
very  local  type,  with  insignificant  value  everywhere  except  in  the 
immediate  area  of  the  inner  and  outer  corners.  The  shear  stresses 
on  the  lower  interface  surface  exhibit  somewhat  different  behavior 
than  the  shear  stresses  on  the  upper  interface  surface.  In  the  case 
of  the  upper  interface,  the  shape  of  the  shear  stress  distribution 
goes  from  a  local  maximum  'negative'  value  at  the  inner  corner  to 
a  local  maximum  'positive'  value  at  the  outer  corner  for  all  cases. 
In  contrast,  at  the  lower  interface  surface,  the  shear  stress 
distribution  for  cases  11  and  10  goes  from  'positive'  to  'negative' 
values  while  cases  7  and  9  go  from  'negative'  to  'positive'  values. 
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Apparently  there  is  some  value  of  the  coefficient  of  thermal 
expansion  of  the  bottom  layer  that  results  in  this  'reversal'  of 
shear  stress  behavior. 


Summary  of  Effect  of  Thermal  Coefficient  of  Expansion  on  System 
Behavior. 

The  study  on  the  effect  of  thermal  coefficient  of  expansion  on 
system  shows  that  overall 

•  The  magnitudes  of  interface  normal  and  shear  stresses  increase 
with  the  the  increase  of  nonuniformity  between  the  thermal 
coefficient  of  expansion  of  the  individual  layers. 

•  Interface  shear  stresses  are  significantly  larger  than 
interface  normal  stresses. 

•  Interface  normal  stresses  may  peeling  or  bearing  depending  on 
the  values  of  the  thermal  coefficients  of  expansion. 

•  Both  the  normal  and  shear  interace  stresses  are  of  an 
extremely  local  nature  and  are  significant  only  in  the  immediate 
area  of  the  corners  of  the  interface. 

It  must  be  pointed  out  that  the  results  obtained  in  the  above  study 
was  for  a  given  geometry  of  the  trimaterial  system.  It  is  expected 
that  these  results  would  hold  qualitatively  for  other  system 
geometries.  The  next  two  study  were  conducted  to  show  the  effect  of 
system  geometrv  on  system  behavior. 
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CHAPTER  4.  THE  EFFECT  OF  MIDLAYER  THICKNESS  ON  SYSTEM  BEHAVIOR 


In  this  chapter,  the  effects  of  nidlayer  thickness  on  system 
behavior  is  investigated.  In  this  study  all  system  properties  (that 
is,  Youngs  modulii  and  thermal  coefficients  of  expansion)  and  all 
system  dimensions  except  midlayer  thickness  are  held  fixed.  A 
series  of  computer  runs  were  conducted  for  different  values  of  the 
midlayer  thickness.  The  following  values  are  fixed.  The  modulus  of 
each  layer  is  100  GPa.  The  coefficients  of  thermal  expansion  of  the 
top,  mid  and  bottom  layers  were  taken  as  100x10'*/ ‘C,  300x10'*/ *C, 
and  10x10'*/ ‘C  respectively.  The  midlayer  thickness  for  the  5  cases 
(5,  12,  13,  14,  and  15)  conducted  are  shown  in  Figure  4.1.  In  the 
5,  12,  13,  14,  15  sequence  the  thickness  of  the  midlayer  decreases 
from  1mm  to  .05mm.  Figure  4.1  also  shows  the  other  dimensions  of 
the  trimaterial  system.  Figures  4.2  through  4.7  show  the  results 
which  are  summarized  below. 

•  For  the  particular  case  of  thermal  coefficients  of  expansion 
used  in  this  study  Figure  4.2  shows  that  the  bearing  normal 
stresses  on  both  the  upper  and  lower  interfaces  decrease  as  the 
thickness  of  the  midlayer  decreases.  Moreover  as  the  thickness 
decreases,  the  distribution  becomes  more  uniform  over  the  interface 
as  well.  In  the  case  of  an  extremely  thin  midlayer  as  might  result 
from  a  thin  adhesive  glue  or  solder,  the  normal  stresses  would  only 
exist  at  the  edges  (corners)  of  the  interface.  Figures  4.2  through 
4.4  present  a  comparison  of  normal  stresses  in  all  of  these  cases 
by  dividing  the  maximum  normal  stress  of  the  particular  case  with 
the  maximum  normal  stress  of  all  the  cases,  that  is,  case  5. 

•  Figures  4.5  shows  the  effect  of  midlayer  thickness  on  the 
interface  shear  stresses.  We  note  that  the  magnitude  of  the 
interface  shear  stress  decrease  as  the  thickness  of  the  midlayer 
decreases.  And  again,  as  the  midlayer  thickness  decreases,  the 
distribution  of  interface  shear  stress  becomes  progressively  more 
and  more  local  in  nature. 

•  Figures  4.6  and  4.7  presenting  comparisons  of  interface  normal 
and  shear  stresses  show  that  the  latter  are  two  to  three  times 
greater.  These  figures  show  that  the  location  of  the  maximum 
interface  shear  stress  on  the  lower  interface  moves  from  the  inner 
edge  to  the  outer  edge  of  the  interface.  On  the  upper  interface 
surface  the  maximum  shear  stress  remained  at  the  inner  edge  of  the 
interface.  The  maximum  normal  stresses  always  occurred  at  the  inner 
edge  for  both  the  upper  and  lower  interface. 

•  Figures  4.2  through  4.7  show  that  all  interface  stresses 
diminish  as  the  thickness  of  the  midlayer  decreases. 

Although  these  results  were  obtained  for  a  particular  fixed  set  of 
thermal  coefficient  properties,  and  layer  lengths,  it  is  expected 
that  these  results  are  valid  qualitatively  for  other  values  of 
thermal  coefficients  of  expansion  and  layer  lengths. 
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KhcW^rC  ap=300xl0'VC  a(.=  lOxlO-^rC  E^  =  Eg=  E(.=  lOOGPa 
t.2  Normal  Stresses  At  The  Interfeces  For  Cases  5,  12,  13,  14,  And  15 


Nomal  Stress/Max  Normal  Stress 


Figure  4.3  2-D  And  3-D  Plots  Of  Normal  Stress/Maximum  Normal  Stress  At  The 
l^jper  Interfece  For  Cases  5, 12, 13, 14,  And  15 
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Normal  Stress/Max  Normal  Stress 


Normal  Stress/Max  Normal  Stress 


Figure  4.4  2-D  And  3-D  Pbts  Of  Normal  Stress/Maximum  Normal  Stress  At  The 
Lower  Inteibce  For  Cases  5, 12, 13, 14,  And  15 
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Position  (mm)  Position  (mm) 

_ Case  5  Case  12  Case  13  Case  14  Case  15 

I  Material  B  Height  Qnm)  1 _ 0;5 _ 0.25 _ OJ _ 0.05 

a^=  100xl0-®/“C  ap=300xl0'^/»C  0^,=  lOxlO^^/^C  E^  =  Eg  =  E(,=  lOOGPa 
Figure4 . 5  Shear  Stresses  At  The  Interfeces  For  Cases  5,  12,  13,  14,  and  15 
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_ Case  5  Case  12  Case  13  Case  14  Case  15 

I  Material  B  Height  ^tim)  1  0.5 _ 0.25 _ OJ _ 0,05 

a^=  100xl0'®/®C  ap=300xl0-^/‘’C  a^,=  10xl0*®/"C  E^  =  Ep=E(.=  lOOGPa 
Figure  4.7  Shear  And  Normal  Stresses  On  The  Midlayer  For  Cases  14  And  15 


CHAPTER  5.  THE  EFFECT  OF  MIDLAYER  LENGTH  ON  SYSTEM  BEHAVIOR 

Here  the  effect  of  nidlayer  length  on  system  behavior  was 
investigated.  The  fixed  values  of  layer  thermal  coefficients  of 
expansion  and  thickness  dimensions  correspond  to  case  5.  That  is, 
the  coefficients  of  thermal  expansion  of  the  top,  mid  and  bottom 
layers  are  100x10  ®/ ‘C,  300x10  */ *C,  and  10x10'*/ *C  respectively,  and 
the  thicknesses  of  each  layer  is  fixed  at  1  mm  each.  The  results  of 
six  cases,  (16,  17,  5,  18,  19,  and  20)  are  compared.  The  sequence  I 

of  these  cases  in  decreasing  magnitude  of  midlayer  length  is  shown 
in  Table  5.1  below.  Figure  5.1  shows  the  geometry  of  these  five 
cases  to  scale. 


1  CASE 

MIDLAYER  LENGTH  mm 

FRACTION  OF  TOTAL  | 

16 

9.0 

• 

VO 

o 

17 

7.5 

.75 

5 

5.0 

.50 

18 

2.5 

.25 

19 

1.0 

.10 

20 

.5 

.05 

Table  5.1  Geometry  for  the  midlayer  length  study 

The  results  of  this  study,  shown  in  figures  5.2  through  5.4,  are  as 
follows. 

•  As  in  all  previous  studies,  the  shear  stresses  are 
significantly  larger  than  the  normal  stresses. 

•  In  all  cases  studied  here,  the  normal  stresses  are  bearing 
stresses. 

•  As  the  midlayer  length  decreases,  the  maximum  interface  normal 
and  shear  stresses  increase  until  the  midlayer  length  becomes  1  mm. 
A  further  decrease  of  the  midlayer  length  decreases  the  stresses 
slightly. 

•  In  all  cases  the  maximum  shear  stress  on  the  upper  interface 
surface  occurs  at  the  inner  corner  point.  The  maximum  shear  stress 
on  the  lower  interface  also  occurs  at  the  inner  corner  point  except 
for  the  least  midlayer  length  case  (case  20)  where  the  maximum 
shear  stress  occurs  at  the  outer  corner  point. 

•  The  maximum  shear  stress  occurs  at  the  inner  corner  point  of 
the  lower  interface  surface  for  the  three  longest  midlayer  length 
cases  (cases  17,  5,  and  18).  For  the  least  midlayer  length  cases 
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Case  20  Lq=  0.5mm 

Figure  5.1  Tri-material  Confgurations  For  Cases  12, 16, 17, 18, 19,  And  20 
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0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9 

NON-DIMENSIONAL  LOCATION  ( x  /  L„ ) 

D 

Figure  5.2  Normal  Stress  At  Inter&ces  As  A  Function  Of  Non-dimensional 
Location  For  Cases  12, 16, 17, 18, 19,  And  20 


a^=  I00xl0'^/“C  ag=300xl0'^/°C  a(,=  lOxlO-^rC  E^  =  Ep=E(,=  lOOGPa 
Figure  5.4  Shear  And  Normal  Stresses  On  The  Midlayer  For  Cases  18,  19,  And  20 


(cases  19  and  20)  ,  the  naxlnuin  shear  stress  occurs  at  the  inner 
corner  point  of  the  upper  interface  surface. 

•  The  naxinuin  normal  stresses  occurs  at  the  outer  corner  points 
of  the  upper  and  lower  interface  surfaces.  The  maximum  normal 
'  normal  stress  on  the  upper  interface  surface  is  10  to  30  percent 

larger  than  the  maximum  normal  stress  on  the  bottom  surface. 

)  Again  it  should  be  pointed  out  that  the  observations  made  above  are 

for  a  particular  trimaterial  system.  The  locations  and  relative 
magnitudes  of  interface  normal  and  shear  stresses  certainly  depend 
on  the  selection  of  the  thermal  coefficients  of  expansion  of  each 
of  the  layers  of  the  trimaterial  system. 
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CHAPTER  6  THE  EFFECT  OF  MONTTMIFORM  TEMPERATURE  FIELD 

This  study  investigated  the  effects  of  a  nonunifom  temperature 
field  on  a  trimaterial  electronic  package  where  the  top  layer  is  a 
silicon  chip,  the  bottom  layer  is  a  substrate  and  a  thin  solder 
midlayer  is  used  to  connect  the  chip  to  the  substrate.  In  a  thesis, 
Sapsai  (4)  undertook  an  exhaustive  study  of  the  thermal  behavior  of 
such  electronic  packages.  A  small  part  of  the  results  obtained  in 
that  study  is  summarized  here. 

The  geometry  and  thermal  properties  of  the  electronic  package  are 
prsented  in  Table  6.1  below. 


Length  (mm) 

Depth  (mm) 

CTE  (/*C) 

E  (GPa) 

CHIP 

5. 

2. 

2.6x10* 

133. 

SOLDER 

.5625 

.2 

29.0X10® 

7.4 

SUBSTRATE 

5.6 

2. 

variable 

variable 

Table  6.1  Geometry  and  Material  Properties  of  Three  Layers 


In  table  6.1  CTE  is  the  coefficient  of  thermal  expansion,  and  E  is 
Youngs  modulus  of  elasticity.  The  Poisson's  ratio's  for  the  Silicon 
chip  and  the  Pb-Sn  solder  are  both  0.3. 

The  particular  materials  comprising  the  electronic  package  is  a 
silicon  chip  in  the  top  layer  with  a  thermal  conductivity 
coefficient  equal  to  150  W/mC  ,,  a  Pb-Sn  solder  in  the  midlayer  with 
a  conductivity  equal  to  63  W/m  .  The  bottom  substrate  layer  varies 
among  four  materials.  Alumina,  Polymide  fiberglass.  Epoxy 
fiberglass,  and  Aluminum  nitride  (AIN).  The  properties  of  these 
four  materials  are  presented  in  Table  6.2  below.  The  Poisson's 
ratio  for  all  three  substrates  is  0.3. 


1  SUBSTRATE 

E  (GPa) 

CONDUCTIVITY 

CTE  (Z'C) 

1  ALUMINA 

262 

18.  W/m'C 

6.0X10  * 

POLYMIDE 

40 

.35  W/m'C 

14.X10® 

EPOXY 

35 

.16  W/m'C 

16.X10  * 

AIN   

339 

230. 

3.3X10* 

The  steady  state  thermal  code  discussed  in  Chapter  2  was  used  for 
the  trimaterial  electronic  package  with  the  geometry  and  material 
properties  given  in  tables  6.1  and  6.2.  The  results  obtained  showed 
that  the  steady  state  temperature  results  was  independent  of  the 
level  of  heat  generation  in  the  chip.  It  was  also  observed  that  the 
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substrate  material  governs  the  temperature  profile  through  the 
package.  The  ratio  of  thermal  conductivities  has  a  significant 
effect  on  system  behavior.  A  characteristic  of  low  thermal 
conductivity  substrates  (i.e.,  kg/k,  <  .0025,  where  subscript  1 
denotes  the  top  chip  layer,  and  subscript  3  denotes  the  bottom 
substrate  layer  is  that  thermal  gradients  occur  in  both  the 
lengthwise  (x)  and  depthwise  (y)  directions.  For  high  thermal 
conductivity  substrates  (i.e.,  .1  <  k,/k3  <  1.6  ),  the  temperature 
field  is  uniform,  that  is,  thermal  gradients  do  not  result.  The 
reasons  for  these  results  are  not  presented  here.  The  interested 
reader  can  refer  to  [4]. 

Figure  6.1  is  a  plot  of  thermal  strain  versus  the  nondimensional 
variable  0y  through  the  thickness  of  the  trimaterial  through  the 
centerline  of  the  three  layers.  The  figure  shows  that  there  is  no 
thermal  strain  gradients  through  either  the  chip  or  solder.  It  also 
shows  the  The  thermal  strain  in  the  solder  is  much  larger  than  the 
thermal  strain  in  the  chip  regardless  of  the  substrate  material,  as 
well  as  much  larger  than  the  thermal  strain  in  the  substrate  in  the 
case  of  the  Aliminum  nitrate  and  the  Alumina  substrates.  In  the 
case  of  the  Epoxy  fiberglass  and  Polymide  fiberglass  substrates, 
the  thermal  strain  increases  monotonically  from  the  bottom  of  the 
substrate  to  the  interface  with  the  solder. 


Figure  6.1  Thermal  strain  through  the  thickness  of  the  trimaterial 
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The  normal  stress  along  the  upper  surface  of  the  solder  midlayer  is 
shovm  in  figure  6.2.  The  figure  shows  that  the  normal  stresses  are 
bearing  in  nature  for  all  substrates.  It  is  also  seen  that  the 
bearing  normal  stresses  on  the  upper  solder  surface  are  much 
larger,  by  a  factor  of  4,  when  Alumina  or  Aluminum  nitrate 
substrates  are  used.  In  all  cases,  the  maximum  normal  stress  occurs 
close  to  the  outer  edge  of  the  solder  layer. 
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Figure  6.2  Normal  Stress  Along  the  Upper  Solder  Surface 

Figures  6.3  and  6.4  show  the  shear  stresses  along  the  upper  and 
lower  surfaces  of  the  solder  midlayer.  These  figures  show  that  on 
the  upper  surface.  Alumina  and  Aluminum  nitride  substrates  achieve 
their  largest  magnitude  at  the  inside  edge,  while  the  shear 
stresses  for  the  Epoxy  fiberglass  and  Polymide  fiberglass  achieve 
their  maximvun  at  the  outer  edge.  Reversal  of  this  behavior  occurs 
along  the  lower  surface.  Shear  stress  reversal  (in  direction) 
occurs  at  various  locations.  In  all  cases  the  shear  stress 
distribution  has  a  constant  slope  along  the  center  two  thirds  span 
of  the  interface.  The  slopes  of  the  shear  stress  distribution 
increase  rapidly  as  the  ends  of  the  interface  are  approached, 
values  at  the  inner  and  outer  edges 
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Figure  6.3  Shear  Stresses  Along  the  Upper  Midlayer  Surface 


Figure  6.4  Shear  Stresses  Along  the  Lower  Midlayer  Surface 
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